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At nonzero density the eigenvalues of the Dirac operator move into the complex plane, while its 
singular values remain real and nonnegative. In QCD-like theories, the singular-value spectrum 
carries information on the diquark (or pionic) condensate. We have constructed low-energy ef- 
fective theories in different density regimes and derived a number of exact results for the Dirac 
singular values, including Banks-Casher-type relations for the diquark (or pionic) condensate, 
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PsJ Smilga-Stern-type relations for the slope of the singular-value density, and Leutwyler-Smilga- 

psj type sum rules for the inverse singular values. We also present a rigorous index theorem for 

non-Hermitian Dirac operators. 
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1. Introduction 

We give a summary of [1], which is very long and sometimes quite technical. Here we present 
the main ideas and most important results. We work in Euclidean space and in the chiral limit. 

2. Dirac eigenvalues and singular values 

In the presence of a chemical potential }i, the massless Dirac operator in QCD and QCD-like 
theories is given by 



with Dy = dy + /Ay, where Ay = A^T„/2 is the gauge field and the T„ are the generators of the 
gauge group in the representation in which the fermions transform. The eigenvalue equation and 
the equation for the singular values are 



where the argument /i is understood if it is not shown explicitly. For /i = we have the trivial 
relation = For nonzero (and real) /i the eigenvalues move into the complex plane, while 
the singular values are always real and nonnegative. Eigenvalues and singular values are unrelated 
for /I / and live on different physical scales. The scale for the eigenvalues is set by the chiral 
condensate at low density [2] and by the BCS gap A at high density [3], while the scale for the 
singular values in QCD-like theories is set by the diquark condensate at all densities, see Sec. 3. 

The operators D(jLt) and D^D have the same zero modes, and the same is true for the operators 
D{—ijl) and DD^ . For a non-Hermitian Dirac operator the index theorem takes the form 



where indD(/i) = dimkerD^ — dimkerD/,. We generically have indD(jLt) = indD(/x)^ unless the 
gauge field is fine-tuned. The proof of (2.3) and further properties of the eigenvalues and singular 
values are given in [1]. 

3. Banks- Casher- type relations for the singular values 

For the derivation of the results in this section it is important that the fermionic measure (in- 
cluding source terms) in the partition function is positive definite, see the discussion in [1, App. C]. 
We restrict ourselves to such cases, which in particular implies an even number Nf of flavors (ex- 
cept for adjoint QCD, where A^^ could also be odd). 

3.1 Two-color QCD 

In two-color QCD the Dirac operator satisfies an anti-unitary symmetry, [CX2K, iD] = with 
C = /7472 the charge conjugation matrix and K the operator of complex conjugation. Adding a 




(2.1) 



D{n)\l/n=KWn 



and D^D(f,„ = 



(2.2) 




(2.3) 
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source term ^jvA-^C75T2/V^ + c.c. with 

/=f ^ -M (3.1) 

to the fermionic Lagrangian and using the Nambu-Gor'kov formaUsm one can show that the parti- 
tion function takes the form 

Z{j) = (det^//2(DtD + /)) = imn+fP"),^^ (3-2) 

\ / YM \ I YM 

where YM stands for Nf = in the average over gauge fields. Taking the logarithmic derivative 
w.r.t. j, followed by the limits four- volume V4 — )• oo and j — )• 0+ (in this order), then yields a relation 
between the scalar diquark condensate and the density Psv(i§ ) of the singular values at the origin, 

{xif^CjsTiIw) = Y^Psv(O) , (3.3) 
see also [4]. Some subtleties of this result and its derivation are discussed in [1]. 

3.2 QCD with isospin chemical potential 

We now consider QCD with Nc > 2 colors, two flavors, and nonzero isospin chemical potential 
jii = Ipt, i.e., we have for the up quark and D{—pL) for the down quark. Adding a source 
term j{dJ^UR — u'^dR) +c.c. to the Lagrangian and using D{pL)^ = — D(— /i) we find 

Z(7) = <det(D^D + /)>Yj^, (3.4) 

which agrees with (3.2) for Nf = 2. After differentiating lnZ(7) w.r.t. j we obtain a relation between 
the pionic condensate and psv, i.e., 

{uysd - dYsu) = 7rpsv(0) . (3.5) 

3.3 Adjoint QCD 

Turning to QCD with adjoint fermions and any number of colors, we add a source term 
^ji/A^CysVA + c.c. and again obtain (3.2). An analogous calculation then yields 

{y^Cy5V) = y^PsM- (3-6) 

4. Low-energy effective theories with diquark sources 

From now on we focus on two-color QCD and distinguish three density regimes that differ in 
their pattern of chiral symmetry breaking and in the number of Nambu-Goldstone (NG) modes. 
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4.1 Low density 

This regime was analyzed in great detail in [5, 6]. Here one starts from the symmetry break- 
ing pattern S\J{2Nf) — Sp{2Nf) at zero density and treats /i and the diquark source as a small 
perturbation. There are Nf{2Nf — 1) — 1 NG modes, parametrized by a field £ = VLdU^ with 
Zrf = diag(/, — /) and U = exp(/7r°r"/2F), where the T" are the generators of the coset space 
SU(2A^y ) / Sp(2A^y ) and F is a low-energy constant (LEC). The leading-order effective Lagrangian 
in the chiral limit is 

772 

ifi^ = — tr(VvIVvE+) - <DLRetr(/E) (4. 1) 




VvL'^ = dvL'^ +H5voCL^B+BL^), (4.2) 

(4.3) 



Jl 

, -J, 



Here, /^/^ are antisymmetric complex matrices of dimension Nf in flavor space, and is another 
LEC equal to the diquark condensate per flavor and handedness at jj. = and without sources, 

<I>L = ^ MWYi) |,-=o,,=o = ^'^) • (4-4) 

For Jr = — Jl = jl there are two types of NG modes [6], 



type 1: mass = J j^h/F^ {Nj-Nf - 1 modes) , (4.5a) 



type 2: mass = ^J J^l/F^ + (2Ai)^ {Nj modes) . (4.5b) 

The type-2 modes become massive for / 0, while the type-1 modes stay massless for j — )• 0. 

4.2 Intermediate density 

At intermediate density jU can no longer be treated as a small perturbation and breaks the 
original SU(2A'^/) symmetry to S\J{Nf)L x S\J{Nf)R x U(1)b. A diquark condensate then breaks 
this symmetry to Sp{Nf)L x Sp{Nf)R. The corresponding NG modes are Ll, Lr G SU(A^/) / Sp{Nf) 
and V G U(1)b, and the total number of NG modes in this regime is Nf{Nf — 1) — 1. The U(1)a 
symmetry is broken explicitly by the anomaly. The effective Lagrangian is 

( Nff^ [\doV\^ + vl\diV\^] + f tr [\doLL\^ + v^\d,LL\^ + {L^R)] 

J -^iRe[Vtr{JLLL-JRLR)] for Nf> 4, (45) 

[ 2/2 [\doV\^ + vl\diV\^] +2^iRe[{jL-jR)V] for Nf = 2, 

where /o,/ and vo,v are LECs corresponding to the decay constants and velocities of the NG 
modes, respectively, and is an LEC similar to (4.4) but with fJ- ^0. All LECs depend on /i. The 
masses of the NG modes are given by 



mA = yJj^l/fl forA = 0,...,Nf{Nf-l)-2, (4.7) 
where = / for A > 1. Note that in this regime all NG modes are massless in the 7 — )■ limit. 
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4.3 High density 



At very high density the U(1)a anomaly is suppressed due to the screening of instantons [7, 8]. 
The U(1)a symmetry of the action is no longer broken explicitly by the anomaly but spontaneously 
by the diquark condensate. Therefore the symmetry-breaking pattern is now [9] 



SU(Af/)L X S\]{Nf)R X U(1)b X U(1)a ^ Sp(^/)l x Sp(iV/)«. 



(4.8) 



The NG modes are the same as at intermediate density, except that there is an additional NG mode 
A G U(1)a which can be considered to be the T]' whose mass has become small. Hence the total 
number of NG modes in this regime is Nf{Nf — 1). The effective Lagrangian is 



H 



Ml 



d^oZ.!" + vg|5,-L|2) + ^ tr (I^oIlP + v2|5,-Il|2) + (L /?) 



-<I>HRetr(7z,LZL - JrRLr) - ^^LtRe {L^Rff^^ 



for Nf > 4, 



[ /2 [\doL\^ + vl\diL\^ + (L o R)] + l^nMhL - JrR) - Mmi^,,Re{L'^R) for Nf = 2, 

(4.9) 



where L = A'^V and R = AV. Similar to Sec. 4.2 we have LECs /o,/ and vq, v as well as Oh, 
all of which depend on /x. The term involving niinst in (4.9) corresponds to the single-instanton 
contribution to the rj' mass, with niinst — )■ as — )• oo, see [1] for a detailed discussion. The masses 
of the NG modes are given by 



type 1: 
type 2: 



(A^^ -Nf-\ modes) , (4. 1 Oa) 

(Imode). (4.10b) 



The type-1 modes stay massless for j — > 0, while the r\' becomes massive as /i is lowered. 

The effective theory at intermediate density can be obtained from the effective theory at low 
density (or high density) by integrating out the NG modes that become massive as the density is 
increased (or decreased). In this way the LECs of the three different regimes can be matched, see 
[1] for details. The domains of validity of the three effective theories, and their overlaps, are also 
discussed in [1]. 



5. Smilga-Stern-type relations for the singular values 

Following the approach of Smilga and Stern for QCD at zero density [10], we can compute the 
slope of the singular- value density at the origin. For technical reasons we start at infinite density. 

5.1 Infinite density 

We now add a more general source term 2VA^C75T2-/V^ + c.c. to the fermionic Lagrangian, 
where J = Jr = —Jl is again an antisymmetric Nf x Nf matrix, which we can decompose as 

N N Y 

J = lY,jAt^=jI + lY.jat" With Ar= A^^.(A^^._l)_l. (5.1) 

A=0 a=l ^ 
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Here, the t^ are the generators of \J{Nf) / Sp{Nf) and the Ja are real parameters with jo 
We can then show the partition function of two-color QCD is given by 

■deti/2(D^D + /7)y,,. = (ndet'^'(^„'+/V); 



Z(7) 



YM 



YM 



For Nf > 4 we define the scalar susceptibility 



KabU) 

for which we obtain after some algebra 



Jim —dj^dj^\nZ{J) 



all ;„=0 



(5.2) 



(5.3) 



KahU) 



lim 



1 

V~4 



£2_ -2 



(^2+^-2)2 



(5.4) 



In the last step we cut off the integral at = A and extracted the part that diverges for j — ^ 0. On the 
low-energy effective theory side we start from (4.9), where at infinite density we can set wijnst = 0. 
Performing a one-loop calculation with a momentum cutoff A we obtain for 7 — )■ 



KabU) ~ 5, 



'ab 



{Nf-4){Nf + 2) 
2Nff^ 

Matching the divergences of (5.4) and (5.5) for j — > 



+ 



yields the slope of Psv at the origin. 



In 



1671^ 



Psv(O) 



{Nf-4){Nf + 2) 



+ ■ 



<T>2 
16712 



(5.5) 



(5.6) 



2Nff^ ■ Nff^r-_ 

This method does not work for Nf = 2, but we expect, based on experience from partially quenched 
chiral perturbation theory [11], that (5.6) is also valid for Nf = 2. 

5.2 Intermediate density 

At intermediate density the calculation on the low-energy effective theory side starts from (4.6) 
and proceeds in a similar way to yield 

'{Nf-4){Nf + 2) 1 



Ps'v(O) 



+ 



167r2 



2Nfr ' Nffif\ 
5.3 Zero density 

At strictly zero density we start from (4.1) with jU = and obtain in a similar way 

{Nf-2){Nf + \) 



(5.7) 



Psv(O) 



5.4 Relation between the three results 



NfF^ 



167r2 ■ 



(5.8) 



At first sight it does not seem possible to interpolate smoothly between the three results for 
Psv(O) at zero, intermediate, and infinite density. This puzzle can be understood by analogy with 
the interpolation between SU(2) and SU(3) chiral perturbation theory, where the strange quark 
mass plays the role of the symmetry -breaking parameter. In our case this parameter is jj.. We find 
(see [1] for details) that at low density the slope is given by (5.7) for E, <^ M^/A and by (5.8) for 
/X^/A <^A, where A ~ F ~ At high density the slope is given by (5.7) for ^ <^ gnr^,/A 
and by (5.6) for gm^,/A <^ ^ <^ gA, where g is the coupling constant. 
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6. Finite-volume analysis 

In a finite volume V4 = L^, tlie three low-energy effective theories constructed in Sec. 4 have 
£-regimes in which the theory becomes zero-dimensional. The corresponding condition is 

— <CL< — , (6.1) 

where rtif is the mass scale of the lightest non-NG particle (see [1, Sec. 5.5] for actual values) and 
mNG is the mass scale of the NG particles, see Eqs. (4.5), (4.7), and (4.10). In these £-regimes 
we have derive Leutwyler-Smilga-type sum rules [12] for the inverse singular values and random 
matrix theories from which microscopic singular- value correlation functions can be derived [1]. 

7. Summary 

We have derived a number of exact results for the singular- value spectrum in QCD-like theo- 
ries without a sign problem. Our results could be used in lattice simulations to obtain the diquark 
condensate at any density, allowing for a numerical test of the conjectured BEC-BCS crossover. 
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